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DYNAMICAL QUARK EFFECTS IN QCD ON THE LATTICE 
RESULTS FROM THE CP-PACS 
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Results of a systematic lattice QCD simulation with two degenerate flavors of sea quarks, identified 
as dynamical u and d quarks, are presented. The simulation was performed on a dedicated parallel 
computer, called CP-PACS, developed at the University of Tsukuba. Clear dynamical quark effects 
are observed in the light hadron mass spectrum and in the light quark masses: In the light hadron 
mass spectrum, major parts of the discrepancy between quenched QCD and experiment are shown to 
be removed by introducing two flavors of dynamical quarks. For the averaged mass of u and d quarks, 
we find m ? (2GeV) = 3.44_ ' 22 MeV using the it and p meson masses as physical input, and for the 

s quark mass, we obtain m^ (2GeV) = 88_ 6 MeV or 90^^ MeV with the K or tf> meson mass as 
additional input. These values are about 20-30% smaller than the previous estimates in the quenched 
approximation. We also discuss the U(l) problem and B meson decay constants. 



1 Introduction 

CP-PACS is a dedicated parallel computer 
designed and developed at the University of 
Tsukuba for simulations in the physics of 
fieldsa. With 2048 node processors intercon- 
nected with a three-dimensional hypercross- 
bar network, the CP-PACS achieves a peak 
performance of 614.4 GFLOPS. Since 1996, 
intensive calculations of lattice QCD have 
been performed on the CP-PACS. Among 
others, the first systematic study including 
both chiral and continuum extrapolations 
was attempted for lattice QCD with two fla- 
vors of dynamical quarks. In this paper, we 
report on the results of these studies, focus- 
ing on the topics of dynamical quark effects 
in QCD. 

We study lattice QCDtl formulated on a 
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4-dimensional hyper-cubic lattice with a fi- 
nite lattice spacing a. Continuum physics is 
defined in the limit of large lattice volume 
and vanishing lattice spacing. Therefore, in 
order to extract predictions for the real world 
from the simulations on finite lattices, we 
have to extrapolate data obtained on a suffi- 
ciently large lattice to vanishing lattice spac- 
ing (the continuum extrapolation). Further- 
more, because the contribution of quarks in 
the calculation is quite computer-time inten- 
sive as we decrease the quark mass, with the 
current computers and current algorithms, 
we also have to extrapolate to the physical 
point of light u and d quarks using data at 
around the s quark mass region {the chiral 
extrapolation). It is important to have good 
control of the systematic errors due to both 
these extrapolations. 

Because of the huge computational power 
required, majority of calculations have been 
made in the quenched approximation, in 
which the effects of dynamical quark loops 
are ignored. As the first project on the CP- 
PACS, we made an extensive simulation of 
quenched QCDH. The quality of extrapola- 
tions and therefore the precision of the final 
hadron spectrum were significantly improved 
over previous studies. From this study, the 
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Figure 1. Quenched light hadron spectrum for 
ground state mesons and baryons in octet and de- 
cuplet representations of flavor SU (3) . 



existence of systematic errors due to the 
quenched approximation was clearly demon- 
strated in the continuum limit. 

Therefore, as the next logical step, we 
then performed a series of "full QCD" sim- 
ulations, in which the effects of dynamical 
quarks are taken into account, on the CP- 
PACSqO. After chiral and continuum extrap- 
olations, clear dynamical quark effects are ob- 
served in the light hadron mass spectrum and 
in the light quark masses. We also found no- 
ticeable effects in B meson decay constants. 

In Sec. 0, we summarize the results for 
the light hadron spectrum from these stud- 
ies. In Sec. 0, light quark masses are dis- 
cussed. Sections |] and [5] are devoted to the 
U(l) problem and B meson decay constants, 
respectively. Conclusions are given in Sec. pi 

2 Light hadron spectrum 

The precise computation of the hadronic 
mass spectrum, directly from the first prin- 
ciples of QCD, is one of the main goals of 
lattice QCD. This provide us with a direct 
and non-perturbative test of the validity of 
QCD as the fundamental theory for strong 
interactions. 

In Fig. [y, the latest results for the light 
hadron spectrum in the quenched approxi- 
mation of QCD are summarized!]. Simula- 
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Figure 2. Continuum extrapolation of vector meson 



m^ and mji- in Nf 



2 and N 



f 



li 



(quenched) QCD, using the K meson mass as input. 



tions are made on four lattices 32 3 x 56 to 
64 3 x 112 with lattice spacings in the range 
a w 0.1-0.05 fm. The spacial lattice size 
was fixed to be about 3 fm, with which the 
finite size effects are estimated to be max- 
imally 0.5% in the spectrum. The u and 
d quarks were treated as degenerate. On 
each lattice, five quark masses, correspond- 
ing to the pseudoscalar-to-vector mass ratio 
mps/my m 0.75-0.4, were studied. The u, d 
quark mass m u d and the lattice spacing a 
were fixed using the experimental values for 
m n and m p as inputs, while the s quark mass 
was fixed either by mjf (if-input) or m$ (</>- 
input). Errors in Fig. H include statistical 
as well as systematic errors from chiral and 
continuum extrapolations, but do not include 
the errors from the quenched approximation. 

From Fig. [I], we see that, although the 
global pattern of the experimental spectrum 
is correctly reproduced, there remain system- 
atic discrepancies of up to about 10% (7 stan- 
dard deviations). The resulting spectrum is 
different depending on the choice of input for 
s quark mass; the K-vcvput or the </>-input. 
These discrepancies and ambiguities are due 
to the quenched approximation. 

Because this limitation of the quenched 
approximation was made clear, the next log- 
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ical step is to perform a "full QCD" calcula- 
tion removing the quenched approximation. 
As the first step towards the realistic QCD, 
we performed a series of QCD simulations 
with two degenerate flavors of sea quarks, 
identified as dynamical u and d quarks, while 
the s quark is treated in the quenched ap- 
proximation (N f = 2 QCD)!!. 

A key ingredient in avoiding a rapid in- 
crease of the computer time is the improve- 
ment of the lattice theory, with which lat- 
tice artifacts are reduced on computationally 
less intensive coarse lattices. We adopted the 
combination of an RG-improved gauge action 
and a "clover" -type improved Wilson quark 
action, and carried out the first systematic 
investigation of full QCD to perform both 
continuum and chiral extrapolations. Our 
preparatory full QCD studyQ shows that this 
action leads already to small lattice artifacts 
at a ~ 0.2 fm. Therefore, we have chosen the 
simulation parameters as summarized in Ta- 
ble [y. The spacial lattice size was fixed to be 
about 2.5 fm for all lattices. 

Recently, we have doubled the statistics 
on the finest lattice at a k O.lfm. All re- 
sults, except for the B meson decay constat- 
nts, presented in this paper are based on this 
full statistics. 

Figure || shows the lattice spacing depen- 
dence of K* and meson masses from Nf = 2 
QCD, compared with the results of quenched 
calculations. For the quenched masses, two 
different data sets are shown: Those denoted 
as "Nf — Standard" are the results of the 



Table 1. Simulation parameters for Ne = 2 QCD on 
the CP-PACS. L s is the spacial size of the lattice. 
On each lattice, four sea quark masses in the range 
mpg/my m 0.8—0.6 were simulated. For each sea 
quark mass, we studied hadrons using five valence 
quarks in the range mpg/mv ~ 0.8-0.5. 



lattice 



a (fm) L s (fm) ^trajectory 



12 3 x 24 


0.215(2) 


2.58(3) 


5000-7000 


16 3 x 32 


0.153(2) 


2.48(3) 


5000-7000 


24 3 x 48 


0.108(2) 


2.58(3) 


4000 



quenched simulation, mentioned before, us- 
ing the standard lattice actioro. Because the 
action used in the full QCD calculation is dif- 
ferent from the original quenched calculation, 
we carried out an additional quenched simu- 
lation using the same improved action as for 
the full QCD runs. The results are denoted 
as "Nf = Improved" in the figure. 

Our data for hadron spectrum confirms 
the expectation that both quenched calcu- 
lations must lead to universal values in the 
continuum limit. The quenched results, how- 
ever, show discrepancies from the experimen- 
tal values, as discussed before. On the other 
hand, when we introduce two flavors of dy- 
namical quarks, the discrepancies are much 
reduced. This means also that the ambigui- 
ties from the choice of input for s quark mass 
are much reduced in Nf = 2 QCD. The re- 
maining small difference might be caused by 
the quenching of the s quark. 

3 Light quark masses 

Although quark massses are the most funda- 
mental parameters of QCD, due to the con- 
finement, it is impossible to measure them 
directly by an experiment. They have to be 
indirectly inferred from hadronic observables 
using a non-perturbative theoretical relation 
between these hadronic quantities and QCD 
parameters. A lattice QCD determination of 
the hadron spectrum provides us with such 
a theoretical relation directly from the first 
principles of QCD. 

Fig. summarizes the lattice spacing de- 
pendence of the average u and d quark mass 
m uc i and the s quark mass m s , in Nf = 2 
full QCD and in quenched QCdI. On the 
lattice, there exist several alternative defi- 
nitions for the quark mass. In the figures, 
they are denoted as VWI (vector Ward iden- 
tity quark masses), AWI (axial- vector Ward 
identity quark masses), etc. See u and for 
details. While different definitions of quark 
masses lead to results that differ at finite lat- 
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Table 2. Light quark masses in the MS scheme at 2 GeV. 

m ud (MeV) m s (MeV) (K-'mput) m s (MeV) (0-input) 



: standard 
: improved 
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4.57±0.18 
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116±3 

notl 

88t 4 6 



144±6 
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Figure 3. Continuum extrapolation of the average u 
and d quark mass m ud and the s quark mass m B in 
the MS scheme at 2 GeV. m s is from the if-input. 
Filled symbols are for Nf = 2 QCD. Quenched results 
with the standard action (qStd) and the improved 
action (qlmp) are shown with thin and thick open 
symbols, respectively. 



tice spacing, they should converge to a uni- 
versal value in the continuum limit. Results 
in Fig. ps clearly demonstrate that this is ac- 
tually the case. 

Values for the light quark masses in the 
continuum limit are summarized in Table |2| 
Errors include our estimates for systematic 
errors from chiral and continuum extrapola- 
tions and renormalization factors. First, we 
note that the two quenched calculations lead 
to universal values, as in the case of the light 
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Figure 4. Continuum extrapolation of the flavor- 
singlet uu + dd meson mass. 



hadron spectrum. However, the quenched 
value for m s differs by about 20% between K- 
input and </>- input. We find that this discrep- 
ancy between the inputs disappears within an 
error of 10% by the inclusion of two flavors 
of sea quarks. 

The most interesting point is that the 
values predicted through Nf = 2 QCD are 
20-30% smaller than those in the quenched 
QCD. In particular, our s quark mass in 
Nf =2 QCD is about 90 MeV, which is sig- 
nificantly smaller than the value w 150 MeV 
often used in hadron phenomenology, and 
almost saturating an estimate of the lower 
bound from QCD sum rules using the pos- 
itivity of spectral functions!!!. On the other 
hand, our result for the u, d to s quark mass 
ratio, Wg /mj^f = 26 ± 2, is consistent with 
24.4 ±^1.5 from one loop chiral perturbation 
theorv 
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4 U(l) problem 



the topological structures. 



The clarification of the mechanism for a large 
rj meson mass is an important issue in QCD. 
Propagators of a flavor non-singlet meson 
consist of a loop of a valence quark propaga- 
tor, while propagators of the flavor singlet rf 
meson have an additional contribution with 
two disconnected valence quark loops. The 
fact that the rf is much heavier than the cor- 
responding non-singlet meson it means that 
the two-loop contribution should exactly can- 
cel the 7r pole of the one-loop contribution, 
leaving the heavy rf pole. This phenomenon 
is considered to be related with the anoma- 
lous violation of the flavor singlet axial U(l) 
symmetry and with the topological structure 
of gauge field configurations. 

The calculation of the two-loop contri- 
bution requires a large amount of computa- 
tions on the lattice. For this reason only lim- 
ited results are available. In an approxima- 
tion ignoring the mixing with the ss state, 
we studied one and two-loop contributions, 
and performed, for the first time, both chi- 
ral and continuum extrapolations!! We ob- 
tain m ufi+d( j = 972 ± 97 MeV for the flavor- 
singlet uu + dd meson. See Fig. [I| In the real 
world, the uu + dd state mixes with the ss 
state to lead to ?7(547) and rf (958) mesons. 
We are extending the study to inspect the 
mixing with the ss state and the relation to 



Table 3. Heavy meson decay constants in MeV. The 
lattice scale was fixed by the p meson mass. Two er- 
rors are statistical and systematic. For B s and D s , 
the s quark mass is fixed from the i^-input; the dif- 
ference between K and </>-input is found to be smaller 
than the systematic error. 







N f = 


N f = 2 


fs d 


Fermilab 


188±3±9 


208±10±11 




NRQCD 


191±5±11 


205±8±15 


/b.. 


Fermilab 


220±2±15 


250±10±13 




NRQCD 


220±5±13 


242±8±17 


fD d 


Fermilab 


218±2±15 


225±14±14 


f Ds 


Fermilab 


250±1±18 


267±13±17 



5 B mesons on the lattice 

The decay constant for the B q meson is de- 
fined by (0|67 M 7 5 (?|B 9 (p)) = if Bq P^ where 
q denotes either d or s quark. The non- 
perturbative determination of fg , and also 
the bag parameters Bb , is quite important 
for a precise determination of CKM matrix 
elements. Therefore, intensive lattice calcu- 
lations have been madetl. 

On the lattice, however, the simulation of 
the heavy b quark is not a trivial extension of 
light quark simulations, because m\, ~ 4 GeV 
is larger than the lattice cutoff ~ 1-4 GeV to 
date. Two methods have been developed to 
simulate heavy quarks on the lattice. One 
is based on a non-relativistic effective theory 
of QCD (NRQCD) defined through an ex- 
pansion in the inverse heavy quark massE9. 
Another employs a relativistic action and 
reinterprets it in terms of a non-relativistic 
Hamiltonian (Fermilab method)E_H. Because 
the both methods include an effective treat- 
ment of heavy quarks, the consistency of the 
results among them should be checked. 

Majority of the lattice studies are done 
in the quenched approximation. Onthc other 
hand, a chiral perturbation theoryEHl suggests 
sizable corrections from dynamical quarks in 
the values of fs ■ The first full QCD cal- 
culations of £n g were made by the MILC 
CollaborationEj using the Fermilab method, 
and by Collins et alH using the NRQCD 
method. In these studies, configurations were 
generated using the staggered sea quarks, 
which is different from the valence light quark 
(the Wilson quarkEj or the clover quarkliJ). 

Using the CP-PACS computer, we stud- 
ied heavy meson decay constants applying 
a consistent formulation for sea and valence 
light quarks, the clover quark, and applied 
both the NRQCD method and the Fermi- 
lab methocfj. Our best estimates of heavy 
meson decay constants for Nf = 2 and 
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Nf = (with improved action) are sum- 
marized in Table g. Because the Fermilab 
method is applicable also for the c quark, 
we also computed fo d and fjj a with this 
method. The fact that our fjj t for Nf — 
2 is consistent with the recent experimen- 
tal results, 285±20±40 MeV (ALEPH0) and 
280±19±44 MeV (CLEO0), is quite encour- 
aging. 

From the table, we see that Fermilab and 
NRQCD methods are consistent with each 
other. We also note that Nf = 2 results 
for B mesons are about 10-15% larger than 
the quenched values, while fjj d and Jd 3 are 
less sensitive to Nf. Increase of the B me- 
son decay constants affects the determination 
of several CKM matrix elements through the 
B q — B q mass difference AM,. Our results 
for fs d and fg 3 are consistent with the hy- 
pothesis that the Wolfenstein parameter p is 
positive. 

6 Conclusions 

We performed the first systematic study of 
lattice QCD with two flavors of dynamical 
quarks. We found that dynamical quark 
effects are quite important in the hadron 
physics. The effect is as large as 20-30% in 
the values of light quark mass and about 10- 
15% in B meson decay constants. Both of 
the shifts has significant implications to phe- 
nomenological studies of the standard model. 
It is urgent to evaluate dynamical quark ef- 
fects in other hadronic quantities, such as the 
bag parameters Bg . It is also important to 
study the effects of dynamical s quark. Fur- 
ther intensive studies on the lattice are un- 
der way to clarify the precise structure of the 
standard model. 
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